Spin squeezing inequalities for arbitrary spin 
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We determine the complete set of generalized spin squeezing inequalities, given in terms of the 
collective angular momentum components, for particles with an arbitrary spin. They can be used 
for the experimental detection of entanglement in an ensemble in which the particles cannot be 
individually addressed. We also present a large set of criteria involving collective observables different 
from the angular momentum coordinates. We show that some of the inequalities can be used to 
detect fc-particle entanglement and bound entanglement. 
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With an interest towards fundamental questions in 
quantum physics, as well as applications, larger and 
larger entangled quantum systems have been realized 
with photons, trapped ions and cold atoms [if. Quan- 
tum entanglement can be used as a resource for certain 
quantum information processing tasks fli], and it is also 
necessary for a wide range of interferometric schemes to 
achieve the maximum sensitivity in metrology 0]- Hence, 
the verification of the presence of entanglement is a cru- 
cial but exceedingly challenging task, especially in an en- 
semble of many, say 10 6 — 10 12 , particles. In such systems, 
typically the particles are not accessible individually and 
only collective operators can be measured. A ubiquitous 
entanglement criterion in this context is the spin squeez- 
ing inequality [3| 
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where N is the number of spin-i particles, Ji 
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for / = x,y,z are the collective angular mo- 



mentum components and jj n ^ are the single spin angular 
momentum components acting on the n th particle. If 
a state violates Eq. ([I]), then it is entangled (i.e., not 
fully separable Q). Such spin squeezed states [j| have 
been created in numerous experiments with cold atoms 
and trapped ions P, @|, and can be used, for instance, in 
atomic clocks to achieve a precision higher than the shot 
noise limit @. 

Recently, after several generalized spin squeezing in- 
equalities (SSIs) for the detection of entanglement ap- 
peared in the literature 0-0] and were used experimen- 
tally [lfl |. a complete set of such entanglement condi- 



tions has been presented in Ref. [ll|. However, all of 
the above mentioned conditions are for spin- 1/2 particles 
(qubits), and so far the literature on systems of particles 
with j > i is limited to a small number of conditions, 
specialized for certain quan tum states or particles with 
a low dimension 0, El , 13 1 . At this point the question 



arises: Could one obtain a complete set of inequalities for 
j > |? Such conditions would be very relevant from the 
practical point of view since in most of the experiments 
the physical spin of the particles is larger than ^ and the 
spin-i subsystems are created artificially. Thus, knowing 
the full set of entanglement criteria for j > |, many ex- 
periments for realizing large scale entanglement could be 
technologically less demanding, and fundamentally new 
experiments could also be carried out. The solution is 
not simple: Known methods for detecting entanglement 
for spin-i particles by spin-squeezing cannot straightfor- 
wardly be generalized to higher spins. For example, for 
j > h , Eq. (fTJ) can also be violated without entanglement 
between the spin- j particles [li| • 

In this Letter, we present the complete set of opti- 
mal spin squeezing inequalities for the collective angu- 
lar momentum coordinates for a system of N particles 
with spin j. We also show how existing entanglement 
conditions for spin-i particles can be transformed into 
entanglement conditions for spin-j particles with j ' > | 
(i.e., qudits with a dimension d = 2j + 1). Finally, we 
present a large set of entanglement conditions for qudit 
systems that involve operators different from the angular 
momentum coordinates, and investigate in detail one of 
the conditions. 

Definitions. The basic idea for the qudit case is that 
besides ji , other single-qudit quantities can also be mea- 
sured. Let us consider particles with d internal states, 
afc for k = 1,2, ...,M will denote single-particle opera- 
tors with the property Tr(afcO;) = CSki, where C is a 
constant. As we will show later, the operators can 
be, for instance, the SU(d) generators for a d dimen- 
sional system. Moreover, for obtaining our generalized 
spin squeezing inequalities, we will need the upper bound 
K for the inequality Sfcli( a l™^) 2 ^ K- 

The ./V-qudit collective operators used in our criteria 

will be denoted by Ak — Y^ n a k^- ^ n tnc qubit case, 
the SSIs were developed based on the first and second 
moments and variances of the such collective operators 
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11]. For j > 1/2, we define the modified second moment 



{ad := (4) - <E( a i n) ) 2 > = E (4 n) 4 m) > (2) 



and the modified variance 

(AA k ) 2 := (AA fe ) 2 - <E(4 n) ) 2 > 



(3) 



In the following, the quantities Eq. ([2]) and Eq. ([3]) will 
be used instead of second moments and variances because 
otherwise it is not possible to obtain tight inequalities for 
separable states [131 ] . 

SSIs for qudits. First, we present a general inequality 
from which the entanglement conditions for the different 
operator sets can be obtained. 

Observation 1. — For separable states, i.e., for states 
that can be written as a mixture of product states [H , 

(N 1) E(A^) 2 ^ £(4) > ~N(N l)K (4) 

kel k0 

holds, where each index set / C {1,2,..., M} defines 
one of the 2 M inequalities. Note that 7 = and 
/ = {1,2, ...,M} are among the possibilities. The proof 
can be found in the Appendix. It is remarkable that the 
bound on the right-hand side of Eq. Q is tight, indepen- 
dent of /, and independent of the particular choice of the 
dfc operators except for the value of K. 

Equation (j4|) is the basis for the entanglement condi- 
tions we present in Obs. 2 and 4. 

Observation 2. — Optimal spin squeezing inequalities 
for qudits. For fully separable states of spin-j particles, 
all the following inequalities are fulfilled 

(Jt) + (Jy) + (J!) < Nj(Nj + l), (5a) 

(A J x f + (A J y ) 2 + (AJ Z ) 2 > Nj, (5b) 

(J 2 k ) + {Jf)-N(N-l)j 2 < (^-l)(AJ m ) 2 , ( 5c ) 

(N - 1) [(AJ k ) 2 + (A J,) 2 ] > (J 2 m )-N(N-l)j 2 ,(5d) 

where k, I, m take all possible permutations of x, y, z. Vio- 
lation of any of the inequalities (J5]) implies entanglement. 
The inequalities ([5|) are a full set for large N in the sense 
that it is not possible to add a new entanglement condi- 
tion detecting other states based on (</&) and (J 2 ). 
Proof. — We applied Observation 1 with {etfc} = 
{jx-,3 Ml jz), K = j 2 and used j 2 + j 2 + j 2 = j{j + 1)1 
the inequalities (O are identi- 
For this case, 
That 



\Jx ! Jyi J 

[3 m 



For j = \ 
cal to the optimal SSIs for qubits 
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the completeness has already been shown 
is, for all values of (Jk) and that fulfill Eqs. (0 

there is a corresponding separable state in the large 
N limit. Direct calculation shows that if a separable 
quantum state £> scp ,i = J2 m Pm.Pm ® p ( m ® •■■ ® pm\ 
where p~m are single-qubit pure states, saturates one 



of the inequalities Eqs. ([5]) for j = i, then the state 

0sep,j = E m Pm u ™ ® w ™ ® — ® w m j saturates the 
same inequality of Eqs. ([5]) for spin-j particles. Here, 

ojffl are single-qudit pure-state density matrices such 
that Tr(prre = ^(^m 3l)- F° r instance, if the first 
state is | + then the second one is | + j) x . Thus 



the proof of completeness of Ref. [ll| can be extended to 
prove the completeness of the criteria Eqs. (JSJ) . □ 
Eq. ([5a| is valid for all quantum states. States maxi- 
mally violating Eq. (|5b|) are angular momentum singlets, 
while for Eq. ((5c|) . for even N, they are symmetric Dickc 



states of the form ( 



N/2\ 
N 



3) 



where Vk denotes all different permutations 17]. 

It is also possible to obtain entanglement conditions 
for spin-J particles from criteria for qubit systems. 
Observation 3. — Let us consider an inequality valid for 
7V-qubit separable states of the form 



l=x,y,z) > COnst., 



(6) 



where / is a concave function of its variables. All of the 
generalized SSIs in the literature have this form. Then, 
the entanglement condition Eq. ^ can be transformed 
to a criterion for a system of N spin-j particles by the 
substitution 



(Ji) 



2.; 



(J?) 



4; 



U(J?))- 



(7) 



Proof. — Let us consider product states of N spin-j par- 
ticles of the form Qj = (EinQ^, and define the quan- 
tities rj = (jj n ^)/j- Then, the first and second mo- 



ments can be rewritten as ( J{) = j 

j 2 E 



» 



and (J 2 



(™)„( m ) 
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The only constraint for the physically 

allowed values for r*j is \r^\ < 1 for all j. Hence, for 
an arbitrary function /, 



t>., 



/({ 2j i^Qi }l=x,y,zi { 47T (Jl )§j }l=x.y,z) 



4,r 



= min/({(J/) e 1/2 }l=x,y,z, {( Jf)g 1/2 }l=x,y,z). 
Ql/2 

If / is a concave function of its variables then we have 
the same minimum for separable states. □ 
Using Observation 3, for instance, the standard spin- 
squeezing inequality Eq. ([T]) from Ref. [3| becomes 



(AJ X ) 
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(Jz 



^ > -■ (8) 



Equation (JS) is violated only if there is entanglement be- 
tween the spin-j particles. Because of the second, non- 
negative term on the left-hand side of Eq. (jHJ , for j > ^ 
there are states that violate Eq. ([1]), but do not violate 
Eq. ([5]). Remarkably, it can be proven that Eq. 
strictly stronger than Eq. ([5]) |17j ]. 



3 



The last application of Obs. 1 is the following. 
Observation 4. — For a system of d-dimensional par- 
ticles, we can define collective operators based on the 
SU(d) generators {g k }tLi with M = <P - 1 as G fe = 



^n=i 9k 1 ' ■ The ^SIs f° r Gk have the general form 



» 



(/V-l)£(AG fe ) 2 -]T(G 2 )> 



kei 



-2N(N-1) 



(d-l) 



(9) 



For instance, for the d = 3 case, the SU(d) generators 

can be the Gell-Mann matrices [18| . 

Proof. — We used Observation 1 with C — 2 and K = 

2(i-j)Bi □ 

Observation 4 presents an abundance of inequali- 
ties. Here, we will analyze in detail Eq. © for / = 
{1, 2, M}. Using Y, k 9l = 2(d+l)(l-£)l |3, Eq. © 
for this case can be rewritten as 



d A -\ 



J2(^Gk) 2 >2JV(d-l). 



(10) 



fc=i 



Direct calculation shows 

noise 



Equation (|10j) is maximally violated by many-body 
SU(d) singlets. Such states appear often in statistical 
physics of spin systems and condensed matter physics 
[201 . They are invariant under operations of the type 
U® N which can be exploited in differential magne- 
tometry [2l| . encoding quantum information in decoher- 
ence free subspaces and sending information independent 
from the reference frame direction 22j . 

Noise tolerance of Eq. U0\) . First, we will ask how 
efficiently Eq. (fTO)) can be used for entanglement detec- 
tion. Let us consider SU(d) singlet states (i.e., states 
with (G 2 ) — 0) mixed with white noise as Q no isy = 

(1 Pnoise ) ^singlet ~t~ Pnoise^pT^- 

that such a state is detected as entangled if p n0 i S c -~ -4+1 ■ 
Thus, the noise tolerance in detecting SU(d) singlets is in- 
creasing with d. Note that Eq. (l5b|) detects a noisy state 
as entangled for an analogous situation if p n oise < gxj ■ 

Eq. H10\) detects k -particle entanglement. The criteria 
presented so far detect any type of non-separability. It 
would be important to find similar criteria that detect 
higher forms of entanglement, that is, fc-entanglement. 
This type of strong entanglement, rather than simple 
non-separability, is needed, for instance, to achieve max- 
imal precision in many interferometric tasks [23| . A pure 
state is said to possess fc-entanglement if it cannot be 
written as a tensor product ® n \ip n ) such that each \ijj n ) 
is a state of at most k — 1 qubits. A mixed state is k- 
entangled if it cannot be obtained mixing states that are 
at most k — 1 entangled 24] . Otherwise the state is called 
(fc — l)-producible. 

While Eq. (fT0|) can be maximally violated by two- 
producible states for j 

j 



21|, it is not the case for 



is needed to violate Eq. (flO|) maximally [15|. Thus, the 
amount of violation of Eq. (|T0|) can be used to detect fc- 
entanglement. 

Observation 5. — For two-producible states the follow- 
ing bound holds 



d A -i 

£ 

fe=i 



(AG fc ) 2 



> 



27V(d- 
27V (d - 



2) 
2) 



for even N, 
2 for odd N. 



(11) 



The violation of Eq. (|TT|) signals 3-particle entanglement. 
Note that for large d the bound in Eq. is very close 
to the bound for separable states in Eq. (fTOj) . The proof 
can be found in the Appendix. 

Eq. hlO\) detects bound entanglement. In Ref. (llj . it 
has already been shown the optimal SSIs for the 3 = \ 
case can detect bound entanglement [25[, i.e., entangled 
states with a positive partial transpose (PPT, [26j]), in 
the thermal states of common spin models. We find nu- 
merically that the criterion Eq. (fTOj) detects bound en- 
tanglement in the thermal state of several Hamiltonians, 
such as for example H = ^ fc G\, even for i > h [171] ■ 

Symmetric states. Next, it is important to ask how 
our entanglement criteria behave for symmetric states, 
as such states naturally appear in many systems such as 
Bose-Einstein condensates of two-state atoms. 
Observation 6. — (i) Symmetric states can violate 
Eq. Q for some / only if pj} n ^ 0, where Tl denotes 
the partial transposition 12 61 ] and the average two-qudit 
density matrix is defined as g av2 = N ^_^ J2 m ^ n Qmn- 
(ii) For symmetric states, if ai~ are the SU(d) generators 
gk, Eq. is equivalent to 



^/V(AG fc ) 2 + (G fe ) 2 >0. 
fee/ 



(12) 



> \. For the SU(d) d-particle entangled state 



For this case, Eq. (TT21 is violated for at least one / 
and some choice of the collective operators if and only 
if gf^ 2 ^ 0- F° r the proof, see the Appendix. 

Implementation. The angular momentum coordinates 
Jfe and their variances can be measured in cold atoms by 
coupling the atomic spin to a light field, and then mea- 
suring the light [6| . The collective spin can be rotated by 

magnetic fields. Measuring the operators X^nOfe™^) 2 can 
be realized by rotating the spin by a magnetic field, and 
then measuring the populations of the j z eigenstates. In 
some cold atomic systems, such operators might also be 
measured directly, as in such systems in the Hamiltonian 
a s z (j^) 2 term appears, where s is the photonic pseu- 
dospin [27]. For the SU(d) generators, the Gk operators 
can be measured in a similar manner, however, SU(2) ro- 
tations realized with a magnetic field are not sufficient. 
For larger spins, it is advantageous to choose the gk op- 
erators to be (|*)<I| + |I)<fc|)/V5, i(\k){l\-\l){k\)/V2imd 
|fc)(fc| 28]. The corresponding collective operators can 
all be measured based an SU(2) rotation within a two- 
dimensional subspace and a population measurement of 
at most two quantum states. 
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In summary, we have presented a complete set of gen- 
eralized SSIs for detecting entanglement in an ensem- 
ble of qudits based on knowing only (J k ) and (J|) for 
k = x,y,z. We extended our approach to collective ob- 
servables based on the SU(d) generators. We showed 
that some of the inequalities can be used to detect k- 
entanglement and bound entanglement. Finally, we dis- 
cussed the experimental implementation of the criteria. 

We thank O. Giihne and Z. Kurucz for discussions. 
We thank the ERC StG GEDENTQOPT, the MICINN 
(Project No. FIS2009-12773-C02-02), the Basque Gov- 
ernment (Project No. IT4720-10), and the National Re- 
search Fund of Hungary OTKA (Contract No. K83858). 

Appendix. — Proof of Observation 1. We consider 
product states of the form |$) = ® n \4>n)- For such 

states, we have {^u) = (Ak) 2 — E„( a L"^) 2 - Hence, 

the left-hand side of Eq. © equals -J2 n ( N ~ 

i)E fce .(4 n) > 2 E^((A.) 2 -E„(4 n) ) 2 ) > 

-En(JV " l)E£ii<4 n) > 2 > ~N(N - l)K. We 
used that (A k ) 2 < N^Ja^) 2 [lljj . □ 
Proof of Observation 5. We will find a lower bound 
on the left-hand side of Eq. JTT]) for N = 2. Let us 
consider first antisymmetric states. We will use that 

E k (G 2 k ) = T.M E fc (l ® 9l) + 2E fc (5fc ® 9k). 

Then, we need that_y] fc §k®gk — 2F — 1 1 where F is the 
flip operator 0,[lJ^Hence, T,k( G t) = 4(d+ 1)(1 - f). 
For the nonlinear part, we have that Efc(fffc)e rod = 
2Tr (e |2 ed) ~ 1 [3 El, and using the Cauchy-Schwarz 
inequality for Efc(#fc ® ® 9k), we obtain a bound 
Efc(^fe) 2 — 4 — §■ Here we used that for antisymmetric 
states, for the reduced single-qudit state Tr(£> 2 cd ) < | 
30 . This leads to Eq. (TTTj) for antisymmetric states. 
For symmetric states the bound on the left-hand side 
of Eq. (fTTj) can be obtained similarly and it is larger. 
Finally, since the equation is invariant under the permu- 
tation of qudits, the variances give the same value for g 
as for \{g + FgF) = P a gP a + P s qP s , where P s and P a 
are the projectors to the symmetric and antisymmetric 
subspaces, respectively. Thus, it is sufficient to consider 
mixtures of symmetric and antisymmetric states. The 
bound for the product of such two-qudit states and of 
single-qudit states for the left-hand side of Eq. ([XT]) can be 
obtained using [AfoSl+KSa)]^^ = (Aa)| i +(Aa) 2 62 . 
Because of the concavity of the variance, the bound is the 
same for mixed 2-producible states. □ 
Proof of Observation 6. Equation Q can 
be rewritten as 52 keI N(AA k ) 2 + (A k ) 2 > 

Efcli(^fe) - N ( N ~ which can be re- 

pressed as Y,kei N ® a k) eav2 ~ ( a k ® l) 2 av2 ) > 
J2 k ( a k ® a k ) em2 - K. From Eq. fll]) for I = it follows 

that E fe ( a fc® a fe)eav2 = jv(jv-i) E fc (^fc) ^ K > while 
the equality holds for symmetric states for the SU(d) 
generators gk [H- We also need that a density matrix 
of a two-qudit symmetric state has a positive partial 



transpose if and only if (O (8 O) — (O ® I) 2 > for every 
O [29(. Hence the statement of Observation 6 follows. 
For qubits, we obtain the results of Ref. Q. □ 
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Supplementary Material 

The supplement contains some derivations to help to 
understand the details of the proofs of the main text. 
It summarizes well-known facts about the quantum the- 
ory of angular momentum and that of SU(d) generators. 
More details will be presented elsewhere [SI]. 

Angular momentum operators. Next, we summa- 
rize the fundamental equations for angular momentum 
operators [S2]. For particle with spin-j we have 



(3x+J« + il)=iO' + i)i- 



(SI) 



Since the angular momentum operators have identical 
spectra, it follows from Eq. (|S1|) that we can write 



1 



Tr(^) = -j(j + l)(2j + l). 



(S2) 



Based on Eq. (|S2j) . we get the constant for the orthogo- 
nality relation 

Tr(jUi) =S k ip(j + l)(2j + l). (S3) 

For the sum of the squares of expectation values we 
have 



k—x,y,z 



(S4) 



For j = i, for all pure states the equality holds for 
Eq. fS4}- 
Finally, 

£ (Ui ® 1 + 1 ® Ji) 2 ) < 2j(2j + 1). (S5) 

l=x,y,z 

Hence, using Eq. (IS 1 [) we obtain 

2j(j + 1) + 2 £ (ji ® j,) < 2j(2j + 1). (S6) 

l=x,y ,z 



Thus, we arrive at the inequality 

£ Oi® ji) <i 2 



(S7) 



Local orthogonal observables. Here we summarize 
the results of Ref. [S3] for Local Orthogonal Observables 
(LOOs, [S4]). For a system of dimension d, these are d 2 
observables \f- such that 



Tr(A fe A ; 



'ki- 



rn 



For a quantum state g, LOOs have the following proper- 
ties 



£(A fc ) 2 = dl, 

k=l 

d 2 

£(A fe ) 2 = Tr(g 2 ) < 1. 



(S9) 
(S10) 



Moreover, based on Ref. [S5] we know that 

d 2 

£A fc 0A fc =F, (Sll) 
fc=i 

where F is the flip operator exchanging two qudits. 

SU(d) generators. Next, we will use the results 
known for local orthogonal observables for SU(d) gen- 
erators. For a system of dimension d, there are d 2 — 1 
traceless SU(d) generators <?& with the property 



(S12) 



Thus, from SU(d) generators gu we can obtain LOOS 
using 



1 



A fe = -^gu 
for k = 1, 2, d 2 — 1, and X d 2 



(S13) 



Vd 



After a derivation similar to that of Ref. [S3] , we arrive 



at 



d 2 -i 



d 2 -l 



(S14) 



fc=i 



1 



£( gfc ) 2 -2^Tr(^)-^j <2^1--). (Si:,) 

fe=l 



1 



/ 1 \ 
£ g k ®g k = 2 I F- -1 J 



(S16) 



Based on Eq. (|S16[) . for bipartite symmetric states we 
have 



df-i 



(J2 9k® 9k) =2 ( +1 



fc=i 



dj' 



(S17) 



while for antisymmetric states we have 



^=1 



^ fffc ®5fc> = 2 [ -1 - - 



(S18) 



k=l 



It is important to stress that the inequalities presented 
are valid for all SU(d) generators, not only for Gell-Mann 
matrices. 

Equations for the collective operators based on 
SU(d) generators. 

Here we present some fundamental relations for the 
collective operators Gk- First of all, the length of the 
vector G — {(Gk)} k=1 is maximal for a state of the 
form |\E')' 8lAr . This can be seen as for such states G = Ng 

-* d 2 — 1 

where g = {(<?/c)*}fe = i , arid knowing that for pure states 
\g\ is maximal. 



6 



For the sum of the squares of Gk we obtain 

fc k n k n^m 

d 2 -l / 1\ 

(S19) 

Here we used Eq. (|S14I) and Eq. (|S16I) . Based on 
Eq. (|S19|) and using (F mn ) > —1, we can write 

£((G fe ) 2 )> — (d + l)(d-N). (S20) 

Note that the bound on the right-hand side of Eq. (|S20|) 
cannot be zero if N < d. For N = d, the sum ^2 k ((Gk) 2 } 
is zero for the totally antisymmetric state for which 
(F m n) = — 1 for all m, n. 
Next, we will show that 

O £(AG fc ) 2 =0. (S21) 

k k 

In order to prove that, one has to notice that 
£ fc (AG fe ) 2 = implies £ fc (AG' fe ) 2 = for any set of 
SU(d) generators G' k [S6]. This also implies (AB) 2 = 
for all traceless observables B. For every traceless D 
one can find traceless B\ and By, such that \B\ , By] = 
iD [S7] and hence (ABi) 2 + (AB 2 ) 2 > \(D)\. Hence, 
^ fe ((Gfc) 2 ) = implies (D) = for all traceless observ- 
ables D [SI]. 

As a consequence of Eq. (|S20I) and Eq. (|S21[) . for 
N < d we have X)fe(^^*;) 2 > 0- Hence, for ci-dimensional 
systems states with less than d particles cannot have 
E fe (AG,) 2 = 0. 

Moreover, for symmetric states we have (F mn ) = +1 
for all m, n, and based on Eq. (IS19[) we obtain 

J2((Gk) 2 ) = —(d-l)(d + N), (S22) 

k 

which is the maximal value for ^2 k ((Gk) 2 ) ■ Similarly, for 



symmetric states, 

E«^) 2 > = £<(G fe ) 2 )-<££G?i n) ) 2 )(S23) 

k k k n 

is also maximal. 

Naturally, these statements are also true for the an- 
gular momentum operators for the j = \ case, as these 
operators, apart from a constant factor, are SU(2) gen- 
erators. 

On the other hand, for the angular momentum oper- 
ators for j > A these statements are not true. In par- 
ticular, (X)fc(^fe) 2 ) i s n °t maximal for every symmetric 
state. 
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